Traffic-driven Epidemic Spreading in Finite-size Scale-Free Networks 
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The study of complex networks sheds light on the relation between the structure and function of complex 
systems. One remarkable result is the absence of an epidemic threshold in infinite-size scale-free networks, 
which implies that any infection will perpetually propagate regardless of the spreading rate. The vast majority 
of current theoretical approaches assumes that infections are transmitted as a reaction process from nodes to 
all neighbors. Here we adopt a different perspective and show that the epidemic incidence is shaped by traffic 
flow conditions. Specifically, we consider the scenario in which epidemic pathways are defined and driven by 
flows. Through extensive numerical simulations and theoretical predictions, it is shown that the value of the 
epidemic threshold in scale-free networks depends directly on flow conditions, in particular on the first and 
second moments of the betweenness distribution given a routing protocol. We consider the scenarios in which 
the delivery capability of the nodes is bounded or unbounded. In both cases, the threshold values depend on 
the traffic and decrease as flow increases. Bounded delivery provokes the emergence of congestion, slowing 
down the spreading of the disease and setting a limit for the epidemic incidence. Our results provide a general 
conceptual framework to understand spreading processes on complex networks. 

PACS numbers: 02.50.Ga,89.75.Fb,89.75.Hc 



I. INTRODUCTION 

Scale-free networks 01 [2, Ht] are characterized by the pres- 
ence of hubs, which are responsible for several striking prop- 
erties for the propagation of information, rumors or infections 
0HI@,I3j[1,IM]- Theoretical modeling of how diseases spread 
in complex networks is largely based on the assumption that 
the propagation is driven by reaction processes, in the sense 
that the transmission occurs from every infected through all 
its neighbors at each time step, producing a diffusion of the 
epidemics on the network. However, this approach overlooks 
the notion that the network substrate is a fixed snapshot of all 
the possible connections between nodes, which does not im- 
ply that all nodes are concurrently active. Many networks ob- 
served in nature 01 Q], including those in society, biology and 
technology, have nodes that temporally interact only with a 
subset of its neighbors Hoi [Tltl . For instance, hub proteins do 
not always interact with all their neighbor proteins at the same 
time IU2I1 . just as individuals in a social network |13j| do not 
interact simultaneously with all of their acquaintances. Like- 
wise, Internet connections being utilized at a given time de- 
pends on the specific traffic and routing protocols. Given that 
transport is one of the most common functions of networked 
systems, a proper consideration of this issue will irreparably 
affect how a given dynamical process evolves. 

Our knowledge of the mechanisms involved in disease 
propagation has improved in the last several years lfl4l[l5i[l6ll . 
Recent works have to some extend surmounted the problem of 
link concurrency through agent-based modeling approaches 
fTill or by explicitly implementing data-driven simulation 



'Electronic address: yamir.moreno@gmail.com 



tools lfl7tl , These very detailed computer models have not 
been yet analyzed in a general theoretical framework. Here 
we introduce a theoretical approach to investigate the outcome 
of an epidemic spreading process driven by transport instead 
of reaction events. To this end, we analyze a paradigmatic 
abstraction of epidemic contagion, the so-called Susceptible- 
Infected-Susceptible (SIS) model OH (see Methods), which 
assumes that contagion occurs through the eventual contact or 
transmission between connected partners that are using their 
connections at the time of propagation. This is achieved by 
considering a quantized interaction at each time step. Mathe- 
matically, we set up the model in a flow scenario where con- 
tagion is carried by interaction packets traveling across the 
network. 

We consider two possible scenarios that encompass most 
of real traffic situations: i) unbounded delivery rate, and ii) 
bounded delivery rate, of packets per unit time. We derive the 
equation governing the critical threshold for epidemic spread- 
ing in SF networks, which embeds, as a particular case, pre- 
vious theoretical findings. For unbounded delivery rate, it 
is shown that the epidemic threshold decreases in finite SF 
networks when traffic flow increases. In the bounded case, 
nodes accumulate packets at their queues when traffic flow 
overcomes the maximal delivery rate, i.e. when congestion 
arises. From this moment on, the results show that both the 
epidemic threshold and the infection prevalence are bounded 
due to congestion. 



II. RESULTS AND DISCUSSION 

We first generate two different types of SF networks. On 
one hand, we construct random SF networks where no corre- 
lations are present using the configuration model 0101- On 
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the other hand, we also generate small-world, SF and highly 
clustered networks — all properties found in many real-world 
networks JH 0] such as the Internet — using a class of re- 
cently developed network models lfl9l l20ll . in which nearby 
nodes in a hidden metric space are connected. This metric 
space can represent social, geographical or any other rele- 
vant distance between the nodes of the simulated networks. 
Specifically, in the model we use, nodes are uniformly dis- 
tributed in a one-dimensional circle and assigned an expected 
degree k from a power law distribution P(k) ~ k 1 . Pairs of 
nodes separated by a distance d are linked with a probability 
r(d; k; k') = (1 + d/d c )~ a , where a > 1 and d c ~ kk' (see 
Methods for further details). In most of our simulations, we 
fixed (k) = 3 and a = 2. Once the networks are built up, the 
traffic process is implemented in the following way. At each 
time step, p = XN new packets are created with randomly 
chosen origins and destinations. The routing of information 
is modeled through even a shortest path delivery strategy or a 
greedy algorithm Il20l 12111 . In the latter, we make use of the 
second class of SF networks and a node i forwards a packet to 
node j in its neighborhood, which is the closest node (in the 
hidden metric space) to the final packet destination. 

To account for link concurrency, we consider that two nodes 
do not interact at all times t, but only when they exchange 
at least a packet. Therefore, we assume that the epidemic 
can spread between nodes every time an interaction takes 
place. This situation is reminiscent of disease transmission 
on air transportation networks; if an infected individual did 
not travel between two cities, then regardless of whether or 
not those cities are connected by a direct flight, the epidemic 
will not spread from one place to the other. In this way, al- 
though a node can potentially interact with as many contacts 
as it has and as many times as packets it exchanges with its 
neighbors, the effective interactions are driven by a second dy- 
namics (traffic). The more packets travel through a link, the 
more likely the disease will spread through it. On the other 
hand, once an interaction is at work, the epidemics spreads 
from infected to susceptible nodes with probability (3. For ex- 
ample, if at time t node i is infected and a packet is traveling 
from node i to one of its neighbors node j, then at the next 
time step, node j will be infected with probability 0. There- 
fore, susceptible and infected states are associated with the 
nodes, whereas the transport of packets is the mechanism re- 
sponsible for the propagation of the disease at each time step. 



A. Unbounded Delivery Rate 

In this situation, congestion can not arise in the system. Fig. 
Q]shows the results for the stationary density of infected nodes 
p as a function of f3 and the traffic generation rate A for SF 
networks. 

The traffic level determines the value of both the epidemic 
incidence and the critical thresholds. We observe the emer- 
gence of an epidemic threshold under low traffic conditions. 
This implies that for a fixed value of A, the epidemic dies out 
if the spreading rate is below a certain critical value /3 C (A). 
More intense packet flows yield lower epidemic thresholds. 



The reason for the dependence of the critical spreading rates 
on A is rooted in the effective topological paths induced by 
the flow of packets through the network. At low values of 
A, there are only a few packets traveling throughout the sys- 
tem, so the epidemic simply dies out because many nodes do 
not participate in the interaction via packets exchanges. As 
A grows, more paths appear between communicating nodes, 
thus spreading the infection to a larger portion of the network. 
Therefore, in traffic-driven epidemic processes the infection 
is constrained to propagate only through links that transmit a 
packet, and thus the number of attempts to transmit the infec- 
tion depends on the flow conditions at a local level, namely, 
on the number of active communication channels at each time 
step. As a consequence, the effective network that spreads the 
infection is no longer equivalent to the complete underlying 
topology. Instead, it is a map of the dynamical process associ- 
ated with packet traffic flow. The conclusion is that the disease 
propagation process has two dynamical components: one in- 
trinsic to the disease itself (/3) and the other to the underlying 
traffic dynamics (the flow). 

To theorize about these effects we next formulate the ana- 
lytical expression for the dependence of the epidemic thresh- 
old on the amount of traffic injected into the system, follow- 
ing a mean-field approach akin to the conventional analysis 
of the reaction driven case. Mathematically, the fraction of 
paths traversing a node given a certain routing protocol I22I1 . 
the so-called algorithmic betweenness, i>£j , defines the flow 
pathways. Let us consider the evolution of the relative den- 
sity, Pk(t), of infected nodes with degree k. Following the 
heterogeneous mean-field approximation |4[], the dynamical 
rate equations for the SIS model are 



d t pk{t) 
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(1) 



The first term in Eq. ((TJ is the recovery rate of infected in- 
dividuals (we set henceforth p = 1). The second term takes 
into account the probability that a node with k links belongs 
to the susceptible class, [1 — Pkit)], and gets the infection 
via packets traveling from infected nodes. The latter process 
is proportional to the spreading probability j3, the probability 
(t ) that a packet travels through a link pointing to an infected 
node and the number of packets received by a node of degree 
k. This, in turns, is proportional to the total number of packets 
in the system, ~ XN, and the algorithmic betweenness of the 
node, b^ ]g . Note that the difference with the standard epidemic 
spreading model is given by these factors, as now the number 
of contacts per unit time of a node is not proportional to its 
connectivity but to the number of packets that travel through 
it. Finally, 0(i) takes the form 



0(t) 



E k b^P(k)p k (t) 
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(2) 



Eq. (dJ has been obtained assuming: (i) that the network is 
uncorrelated P{k'\k) = k'P{k')/(k), and (ii) that the algo- 
rithmic flow between the classes of nodes of degree k and k' 
factorizes ~ b^ lg b^ lg . Although no uncorrelated networks 
exist, this approximation allows us to identify the governing 
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parameters of the proposed dynamics. The second approxima- Within our formulation, the classical result |01 



tion is an upper bound to the actual value of the ti^ , whose 
mathematical expression is, in general, unknown. The valid- 
ity of the theory even with these approximations is notable as 
confirmed by the numerical simulations. 

By imposing stationarity [dtPk(t) = 0], Eq. (Q]) yields 



Pk 



(3) 



from which a self-consistent equation for is obtained as 

1 ^ (b^) 2 P(k)f3\NG 



e = 



(4) 



The value 6 = is always a solution. In order to have a 
non-zero solution, the condition 



1 
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>1 (5) 
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must be fulfilled, from which the epidemic threshold is ob- 
tained as 
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below which the epidemic dies out, and above which there is 
an endemic state. In Fig.|2]a comparison between the theoret- 
ical prediction and numerical observations is presented. Here, 
we have explicitly calculated the algorithmic betweenness for 
the greedy routing as it only coincides with the topological 
betweenness for shortest paths routing. The obtained curve 
separates two regions: an absorbing phase in which the epi- 
demic disappears, and an active phase where the infection is 
endemic. 

Equation ^ is notably simple but has profound implica- 
tions: the epidemic threshold decreases with traffic and even- 
tually vanishes in the limit of very large traffic flow in fi- 
nite systems, in contrast to the expected result of a finite- 
size reminiscent threshold in the classical reactive-diffusive 
framework. Admittedly, this is a new feature with respect to 
previous results on epidemic spreading in SF networks. It is 
rooted in the increase of the effective epidemic spreading rate 
due to the flow of packets. This is a genuine effect of traffic- 
driven epidemic processes and generalizes the hypothesis put 
forward in the framework of a reaction-diffusion process [15] 
on SF networks. It implies that an epidemic will pervade the 
(finite) network whatever the spreading rate is if the load on it 
is high enough. Moreover, Eq. © reveals a new dependence. 
The critical threshold depends on the topological features of 
the graph, but at variance with the standard case, through the 
first two moments of the algorithmic betweenness distribu- 
tion. As noted above, the algorithmic betweenness of a node 
is given by the number of packets traversing that node given 
a routing protocol. In other words, it has two components: a 
topological one which is given by the degree of the node and 
a dynamical component defined by the routing protocol. 
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can be obtained for a particular protocol and traffic conditions, 
although we note that the microscopic dynamics of our model 
is different from the classical SIS. To see this, assume a ran- 
dom protocol. If packets of information are represented as w 
random walkers traveling in a network with average degree 
(k), then under the assumption that the packets are not in- 
teracting, it follows that the average number of walkers at a 
node i in the stationary regime (the algorithmic betweenness) 
is given by J231I3 K\z ~ TJJkJ w - The effective critical value 
is then (0X) C =< k > 2 /(< k 2 > w), that recovers, when 
u> = (k), the result in Eq. (0. 

Our results are robust for other network models and differ- 
ent routing algorithms. We have also made numerical simula- 
tions of the traffic-driven epidemic process on top of Barabasi- 
Albert and random SF networks implementing a shortest paths 
delivery scheme. In this case, packets are diverted follow- 
ing the shortest path (in the actual topological space) from the 
packets' origins to their destinations. The rest of model pa- 
rameters and rules for epidemic spreading remain the same. 
Figure [3] shows the results obtained for random SF networks. 
As can be seen, the phenomenology is the same: the epidemic 
threshold depends on the amount of traffic in the network such 
that the higher the flow is, the smaller the epidemic threshold 
separating the absorbing and active phases. On the other hand, 
for processes in which the delivery of packets follows a short- 
est path algorithm, Eq. (|6]l looks like 
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where fe top is the topological betweenness. 



B. Bounded Delivery Rate 

Equation (H)) allows us to investigate also the equivalent sce- 
nario in the presence of congestion. Let us consider the same 
traffic process above but with nodes having queues that can 
store as many packets as needed but can deliver, on average, 
only a finite number of them at each time step. It is known 
that there is a critical value of A above which the system starts 
to congestJH] 



A, 



(iV-1) 



(9) 
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Equation (0 gives the traffic threshold that defines the onset 
of congestion, which is governed by the node with maximum 
algorithmic betweenness b'* lg . Substituting (0 in (O we obtain 
a critical threshold for an epidemic spreading process bounded 
by congestion. Increasing the traffic above A c will gradually 
congest all the nodes in the network up to a limit in which the 
traffic is stationary and the lengths of queues grow without 
limit. 
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To illustrate this point, let us assume that the capacities 
for processing and delivering information are heterogeneously 
distributed lEM l26l l27ll so that the larger the number of paths 
traversing a node, the larger its capability to deliver the pack- 
ets. Specifically, each node i of the network delivers at each 
time step a maximum of [~Cj = 1 + k%] packets, where 77 is a 
parameter of the model. In this case, the critical value of A in 
Eq.© is multiplied by the maximum delivery capacity l25ll . 
Moreover, without loss of generality, we will explore the be- 
havior of the model in random SF networks where the routing 
is implemented by shortest paths b a i g = b top ~ k v , being v 
usually between 1.1 and 1.3 JUt]. The previous assumption 
for the delivery capability thus allows to explore as a func- 
tion of 77 the situations in which the delivery rate is smaller 
or larger than the arrival rate (defined by the algorithmic be- 
tweenness). Phenomenologically, these two scenarios corre- 
spond to the cases in which the traffic is in a free flow regime 
(if 77 > v) or when the network will congest (if 77 < v). We 
also note that the adopted approach is equivalent to assume a 
finite length for the queues at the nodes. 

Figure [4] shows the results for the evolution of the average 
density of infected nodes, p, as a function of the spreading rate 
/3 and the rate at which packets are generated A for two differ- 
ent values of 77 using a shortest path delivery scheme on top of 
random SF networks. For 77 = 0.8, the epidemic incidence is 
significantly small for all values of the parameters A and (3 as 
compared with the results obtained when the rate of packets 
delivery is unbounded. On the contrary, when 77 = 1.7 the 
phase diagram is qualitatively the same as for the unbounded 
case, including the result that the epidemic incidence vanishes 
when A is large enough. A closer look at the dynamical evo- 
lution unveils an interesting, previously unreported, feature — 
when the rate at which packets are delivered is smaller than the 
rate at which they arrive, the average value of infected nodes 
saturates beyond a certain value of the traffic flow rate A. This 
effect is due to the emergence of traffic congestion. When the 
flow of packets into the system is such that nodes are not able 
to deliver at least as many packets as they receive, their queues 
start growing and packets pile up. This in turns implies that 
the spreading of the disease becomes less efficient, or in other 
words, the spreading process slows down. The consequence 
is that no matter whether more packets are injected into the 
system, the average level of packets able to move from nodes 
to nodes throughout the network is roughly constant and so is 
the average level of infected individuals. 

Figure|5]illustrates the phenomenological picture described 
above. It shows the epidemic incidence p for a fixed value of 
P = 0.15 as a function of A for different values of 77. The 
figure clearly evidences that congestion is the ultimate reason 
of the behavior described above. Therefore, the conclusion 
is that in systems where a traffic process with finite delivery 
capacity is coupled to the spreading of the disease the epi- 
demic incidence is bounded. This is good news as most of the 
spreading processes in real-world networks involves different 
traffic flow conditions. Further evidence of this phenomenol- 
ogy is given in Fig. [6] where we have depicted the epidemic 
threshold as a function of A for two different values of 77, less 
and greater than v. When 77 < v congestion arises, and the 



contrary holds for 77 > v where the diagram is equivalent to 
that of unbounded traffic. The onset of congestion determines 
the value of j3 above which congestion starts. It is clearly vi- 
sualized as the point beyond which the power law dependence 
in Eq. © breaks down. The plateau of j3 c corresponds to the 
stationary situation of global congestion. 



m. CONCLUSIONS 

In our study, we have computed both analytically and nu- 
merically the conditions for the emergence of an epidemic 
outbreak in SF networks when disease contagion is driven 
by traffic or interaction flow. We demonstrate that epidemic 
thresholds are determined by contact flows and that the epi- 
demic incidence is strongly correlated to the emergence of 
epidemic propagation pathways that are defined and driven 
by traffic flow. Two general phenomenological scenarios are 
identified. When the rate at which packets are delivered is 
large enough or it is not bounded, the epidemic threshold 
is inversely proportional to the traffic flow in the system. 
This simple law, which also encompasses a dependency with 
the topological properties of the graph as given by the algo- 
rithmic betweenness has remarkable consequences. One of 
them is that the epidemic threshold vanishes in the limit of 
high traffic conditions, even for finite SF networks and for 
any routing algorithm. Moreover, within this scenario, the 
new approach recovers previous results as a particular case of 
reaction-diffusion spreading. New phenomenological proper- 
ties appear, in the unified framework, when the system has a 
finite capability to deliver packets and congestion is possible. 
In this case, what determines the epidemic incidence is the 
capability of the network to avoid congestion. If congestion 
arises, the number of contacts between the system elements 
decreases, leading to a less efficient spreading of the disease 
and therefore to a significant reduction of the average num- 
ber of infected individuals. Additionally, the value of the epi- 
demic threshold depends on the delivery rate — the smaller the 
delivery rate, the larger the epidemic threshold. Interestingly, 
our study points out that congestion appears to have a posi- 
tive effect reducing the critical spreading rate for an epidemic 
outbreak to occur and its size. 

Finally, we point out that the new approach presented here 
constitutes an interesting framework to address related prob- 
lems. For instance, in the context of air- transportation net- 
works |[l7l l29ll . traffic-driven mechanisms of disease propa- 
gation could explain the observed differences in the impact 
of a disease during a year l3(ill . One might even expect that, 
due to seasonal fluctuations in flows, the same disease could 
not provoke a system-wide outbreak if the flow were not high 
enough during the initial states of the disease contagion. The 
same reasoning implies that air traffic flow restrictions is an 
efficient way to content a pandemic spreading. Incorporat- 
ing the traffic-driven character of the spreading process into 
current models has profound consequences for the way the 
system functions. Also the theory could help designing new 
immunization algorithms or robust protocols; one in particu- 
lar being quarantining highly sensitive traffic nodes or equiv- 
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alently control the onset of congestion in networks (control 
of air traffic flow). Obviously our approximation lacks of the 
level of detail required to assess more realistic epidemiologi- 
cal models, however it provides a simple framework to antic- 
ipate most of the general features to be observed in detailed 
computational agent-based models. 



it through node j that is the closest to I in the hidden met- 
ric space. Alternatively, one can implement a shortest path 
algorithm in which packets are delivered through nodes that 
are closer to their destinations in the network (i.e., not in the 
hidden metric space). Our results are insensitive to the two 
routing protocols implemented. 



IV. METHODS 
A. Structure 

Scale-free networks are generated by assigning all nodes 
a random polar angle 8 distributed uniformly in the interval 
[0, 2tt). The expected degrees of the nodes are then drawn 
from some distribution x(k) and the network is completed 
by connecting two nodes with hidden coordinates (8, k) and 

(0',fc') with probability r(6,k,6',k') = (l + ^jff)^, 

where rf = (a — l)/2(fc), d(8,8') is the geodesic distance 
between the two nodes on the circle, and (k) is the aver- 
age degree. Finally, choosing x(k) = (7 — l)fcQ _1 fc _7 , 
k > ko = (7 — 2)(k)/(j — 1) generates random networks 
with a power law distribution with exponent 7 > 2 . 

B. Routing Dynamics 

At each time step, p = XN new packets are created with 
randomly chosen origins and destinations. For the sake of 
simplicity, packets are considered non-interacting so that no 
queues are used. The greedy routing delivers packets such 
that node i, holding a packet whose destination is I, dispatches 



C. Epidemic Dynamics 

We have implemented the Susceptible-Infected-Susceptible 
model in which each node can be in two possible states: 
healthy (S) or infected (I). Starting from an initial fraction 
of infected individuals po = Io/N, the infection spreads in 
the system through packet exchanges. A susceptible node has 
a probability (3 of becoming infected every time it receives a 
packet from an infected neighbor. We also assume that in- 
fected nodes are recovered at a rate /1, which we fix to 1 for 
most of the simulations. After a transient time, we compute 
the average density of infected individuals, p, which is the 
prevalence of disease in the system. 
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FIG. 1 : Dependence of epidemic incidence on traffic conditions for unbounded delivery rate. The density of infected nodes, p, is shown as a 
function of the spreading rate f3 and the intensity of flow A in SF networks. Flow conditions (controlled by A) determine both the prevalence 
level and the values of the epidemic thresholds. Increasing the number of packets traveling through the system has a malicious effect: the 
epidemic threshold decreases as the flow increases. Each curve is an average of 100 simulations starting from an initial density of infected 
nodes po = 0.05. The results corresponds to the greedy routing scheme and the network is made up of 10 3 nodes using the model in l20ll . The 
remaining parameters are a = 2, 7 = 2.6 and (k) = 3. 
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FIG. 2: Comparison between numerical and theoretical critical points. Log-log plot of the critical thresholds, /3 C , as a function of the rate at 
which packets are injected into the system, A. Two regions are differentiated: an active and an absorbing phase as indicated. The solid line 

corresponds to Eq. l[6} with jtttt? = 0.154. The agreement is remarkable even though Eq. ® is derived using a mean field approach. The 

\ alg' 

underlying network, infection spreading mechanism and routing protocol are the same as in Fig.[T] Each curve is an average of 10 2 simulations. 
Remaining parameters are the same as in Fig.Q] 
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FIG. 3: Density of infected nodes, p, as a function of traffic flow (determined by A) and the epidemic spreading rate (3 for random SF networks 
and a shortest paths routing scheme for packets delivery. Each point is the result of 100 averages over different networks and initial conditions. 
The network has a degree distribution with an exponent 7 = 2.7. 




FIG. 4: Dependence of epidemic incidence on traffic conditions for the bounded delivery rate scenario. The figures show the density of 
infected nodes, p, as a function of the spreading rate (3 and the intensity of flow A in random SF networks. The top panel corresponds to the 
case 77 = 0.8, for which the epidemic threshold is determined by congestion. When the whole network gets congested, no matter whether the 
value of A is further increased, j3 c remains constant. In addition, the epidemic incidence saturates. The bottom panel shows the results for 
rj — 1.7. In this case, the epidemic threshold vanishes when A grows and Eq. l[6} holds, thus resembling the case of unbounded delivery rate 
(Fig.[Ul. Each curve is an average of 100 simulations starting from an initial density of infected nodes po = 0.05. The network is a random SF 
network made up of 10 3 with 7 = 2.7 and (k) — 3. 
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FIG. 5: Epidemic incidence in traffic-driven epidemic processes with bounded delivery rate. The figure represents the average fraction of 
infected nodes p as a function of A for different delivery rates at fixed j3 = 0.15. The curves depart from each other when congestion arises 
and the epidemic incidence saturates soon afterward. The rest of parameters are those of Fig.[4j 
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FIG. 6: Epidemic thresholds as a function of A for two values of r\. The onset of congestion (the arrow in the figure) marks the point, 
A c ~ 0.150, at which the curve for r\ = 0.8 departs from Eq. i.e., when the power law dependence breaks down. Soon afterwards 
congestion extends to the whole network leading to a bounded (from below) epidemic threshold. On the contrary, when the delivery rate 
is large enough (as in the case of r\ = 1.7), Eq. ijfj} holds for all values of A, thus resembling the unbounded delivery rate case. Network 
parameters are those of Fig. [4] 



